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Abstract 

For the reference triangle or tetrahedron T, we study the stability properties of the L 2 (T)-projection 
IIjv onto the space of polynomials of degree N. We show ||IIjv"it||^2 (^7-) — C|l n llL 2 (T) IMI-fPCn - This 
implies optimal convergence rates for the approximation error \\u — I1jv"'«||l 2 ( ( 97 _ ) f° r an u ^ H k (T), 
k > 1/2. 

1 Introduction and main results 

The study of polynomials and their properties as the polynomial degree tends to infinity has a very 
long history in numerical mathematics. Concerning approximation and stability properties of various 
high order approximation operators, the univariate case is reasonably well understood (in the way of 
examples, we mention the monographs J2BJ for orthogonal polynomials and [U] for issues concerning 
approximation). One-dimensional results can often be generalized to the case of the <i-dimensional 
hyper cube by tensor product arguments. Indeed, a significant number of results is available for this 
case, and we point the reader to the area of spectral methods [TJ|3] and to [15] and references therein. 
The situation is less developed for simplices, possibly due to a presumed lack of product structure. It is 
the purpose of this note to contribute to this field by studying the stability properties of the polynomial 
L 2 -projection on triangles and tetrahedra, paying special attention to its trace on the boundary. Our 
main result is Theorem 11.11 below on the stability of the boundary trace of the L 2 -projection. It 
generalizes known results for tensor product domains to the case of triangles/tetrahedra: Theorem ll.il 
is the analog of Lemma 4.2] (and correspondingly, Corollary 11.21 is the analog of [144 Rem. 4.3] 
and [17} Lemma 3.5]). We mention that the two-dimensional case of triangles of Theorem 11.11 was 
proved in [22]. Independently, closely related results have recently been obtained in [5J and [T2J. 
The novelty of the present work over [6] is that, in the language of Corollary 11.21 below, we extend 
the approximation result of [B] from s > 1 to s > 1/2. The elegant proof of Corollary 11.31 given 
in [121 Lemma 6.4] relies heavily on Corollary 11.21 the analysis presented here provides the tools for 
an alternative proof of Corollary 11.31 which we outline as well. 

The stability result of Theorem 1 1 . 1 1 has applications in the analysis of the hp- version of discontinuous 
Galerkin methods (/ip-DGFEM) as demonstrated in [23]. More generally, simplicial elements are, 
due to their greater geometric flexibility as compared to tensor product elements, commonly used in 
high order finite element codes so that an understanding of stability and approximation properties of 
polynomial operators defined on simplices could be useful in other applications of high order finite 
element methods (/ip-FEM) as well. We refer to [7,8,18,22,23] for various aspects of /ip-FEM. 
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To fix the notation, we introduce the reference triangle T 2 , the reference tetrahedron 7~ 3 as well 



as the reference cube S d by 

T 2 := {(x,y) G M 2 : -1< x < 1 , -Ky<-x}, (1.1a) 

T 3 := £l 3 :-l<i,|/,Kl, x + y + z < -1}, (1.1b) 

S d :=(-l,l) d , d€ {1,2,3}. (1.1c) 



Throughout, we will denote by Vn the space of polynomials of (total) degree N. We then have: 

Theorem 1.1. Let T be the reference triangle or tetrahedron and denote by II jv : L 2 (T) — > Vn the 
L 2 (T) -projection onto the space of polynomials of degree N. Then there exists a constant C > 
independent of N such that 

\\ u Nu\\ 2 L 2 (dT) < C \\u\\tf(j) ||w||hi(T) Vu g hX (T). (1.2) 

In particular, therefore, 

\\H-Nu\\ L 2 {dT) <C\\u\\ B i/2 {T) V«G Bl /2 (T), (1.3) 

where the Besov space B^ 2 {T) = (L 2 (T), H l {T))\/2,i is defined as an interpolation space using the 
so-called real method (see, e.g., for details). 

Corollary 1.2. Let T be the reference triangle or tetrahedron. Then for every s > 1/2 there exists a 
constant C s > such that 

\\u-UNu\\^ m ^CsiN + l^-^WuWHsrr) \/ugH s (T). (1.4) 

The following result, which generalizes the analogous result on hyper cubes of Thm. 2.2], is 
derived in |12[ Lemma 6.4] with an integration by parts argument, a polynomial inverse estimate, and 
Corollary II. 21 with s = 1. 

Corollary 1.3 ( [12j Lemma 6.4]). LetT be the reference triangle or tetrahedron. Then there exists 
a constant C > such that for all N € N 

||n^|| H1(r) < C(N + l) l ' 2 \\u\\ m{T) Vu G H\T). (1.5) 

We will only prove Theorem 11.11 and Corollary 11.21 for the 3D case at the end of Section [SJ We refer 
to the Bachelor Thesis |29j for the 2D version. We will present the key steps for an alternative proof 
of Corollary 11.31 at the end of Section [6l 

1 /2 _____ 

Remark 1.4. The regularity requirement u G B 21 (T) in the stability estimate f j 1 . 3 j) is essentially 
the minimal one, if uniformity in N is sought. To see this, we first note that taking the limit N — > oo 
reproduces the known result (see, e.g., [____, Chap. 32]) ||7o^||l 2 (T) — C\\ u \\ R i/2 m , where 70 is the trace 

operator. Next, we recall that for any e > we have iJ 1 / 2+e (T) C B\ /2 (T) C # 1/2 ~ e (T) but the 
trace operator 70 cannot be extended to a bounded linear operator H l l 2 ~ e {T) — > L 2 (dT). ■ 

1.1 Outline of the proof of Theorem 11.11 

In the tensor-product case of squares and hexahedra, the arguments leading to Theorem 11.11 and 
Corollary 11.31 can be reduced to a one-dimensional setting and have been worked out in [T__] — we 
recapitulate the key steps later in this subsection. Most of this reduction to one-dimensional settings 
is also possible for simplices, where Jacobi polynomials come into play instead of the simpler Legendre 
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or Gegenbauer/ultraspherical polynomials. Let us highlight some reasons why a reduction to one- 
dimensional situations is possible and what the key ideas of the proof of Theorem 11.11 are: 

(I) The first basic tool for a dimension reduction is a classical one, the so-called Duffy transfor- 
mation (see (|4.1|) below), which maps the simplex into a hyper cube. As noted already by |!U|, ll8| , [Tn] 
orthogonal polynomials on the simplex (see Lemma l4.3p can be defined in the transformed variables 
through products of univariate Jacobi polynomials, which expresses the desired reduction to one- 
dimensional settings. These Jacobi polynomials arise since the transformation from the simplex to 
the cube changes the classical (unweighted) Lebesgue measure to a weighted measure. The situation is 
technically more complicated than the tensor-product setting: first, Jacobi polynomials arise (due to 
the weights) in contrast to the more common Gegenbauer/ultraspherical polynomials for the tensor- 
product setting second, the weight in the Jacobi polynomials is not fixed so that the dependence 
on both the polynomial degree and the weight needs to be tracked explicitly (cf. the definition of 
the orthogonal polynomial on the simplex in Lemma |4,3|) . Nevertheless, the Jacobi polynomials are 
classical orthogonal polynomials, and one can draw on a plethora of known properties for the purpose 
of both analysis and design of algorithms. We mention in passing that these observations have been 
made and exploited previously in different connections, for example, in the works [2"|lo]ll81l2"C)]. 

(II) A second ingredient to the reduction to one-dimensional problems arises from the fact that 
we aim at trace estimates in Theorem ll.il Such estimates hark back to techniques associated with the 
names of Gagliardo and Nirenberg and are closely connected to the ID Sobolev embedding H 1 C L°°. 
For example, for the half space := {(x',Xd) \x' G R^ -1 ,^ > 0}, the trace estimate can be cast 

in multiplicative form as 0)||^ 2 , Kd _ 1 ^ < [|w[|x,2(r<») IMIz 2 (R d ) + ll^o"^^) where one has some 

freedom to choose the vector (field) a as long as it is not tangential to the plane M. d ~ l x {0} (the 
classical proofs usually take a to be a constant normal to x {0}). This freedom to suitably 

choose the vector field a is exploited in the proof of Corollary 11.21 given in [6] , where a is a vector 
field pointing from one face of the reference tetrahedron to the vertex opposite. Our analysis here and 
that of [29] effects a similar thing: It performs the analysis on the reference cube 5 3 for the Duffy- 
transformed function u and singles out d m u, which is closely related to the directional derivative d a u 
on the simplex. 

(III) Given that we study the L 2 -projection, the observation | (II) shows that we need a connection 
between the expansion of u in orthogonal polynomials on the triangle/tetrahedron and the expansion 
of d a u. This is at the heart of the analysis of [6]. Likewise it is the key step in [29J and the present 
article. Working in the transformed variables on the cube S d and denoting the transformed function 
u, we have to relate the expansion of u to that of d V3 u. Ultimately, the issue is to understand the 
relation between the sequences (u q )'^ =0 and (bg)^L , where the terms u q and b q are the coefficients of 

the function u = u q P q a ' ^ and its derivative u' = Y^ q b q P q a '°^ in terms of the expansion in series of 

Jacobi polynomials P q a '°\ A technical complication over the tensor product case is that a family of 
expansions (parametrized by a) has to be considered and that the dependence on a has to be tracked. 
It is easy to express the coefficients u q in terms of the coefficients b q — see Lemma 13.31 for the general 
case and (|1.6|) for the special case a = 0. It is harder to control the coefficients b q in terms of the 
coefficients u q . In the present work, this is achieved in Lemma 13.61 and discussed in more detail in the 
following point |(IV)[ 

(IV) The refinement of the present analysis and [29] over [6] is the multiplicative structure of 
the estimate. This results from a refined connection between the expansion of a function and its 
derivative. Let us review the origin of the multiplicate estimate (|1.2p for tensor product domains 
as given in [14, Lemma 4.1]), since similar ideas underlie our arguments. For the case N > 2 and 
T 1 := (—1,1), we aim to show |(ILv , u)(l)| 2 ^5 ll n llL 2 (T 1 )ll' u ll-ff 1 (T 1 )- By the ID Sobolev embedding 
theorem already mentioned above we have ||«|| £00(7-1) ^5 II u IIl 2 (T 1 ) ll^ll/J" 1 (7 -1 ) • ^ therefore suffices to 
establish the inequality for u — H^u. Following [171 Lemma 3.5], we expand u and its derivative u' in 
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Legendre series (we write, as is common, L q = P q ' ^): 



q=0 q=0 

with coefficients u q and b q explicitly given by u q = j-j-i u(x)L q (x) dx and b q = J^-i u'(x)L q (x) dx. 
Orthogonality properties of the Legendre polynomials imply (see Lemma 13.31 below with a = 0) 

u _Azl h±L_ a>2 (16) 

Uq ~2q-l 2g + 3' 9 ~ 1 bj 



Since L q (l) = 1 for all q G No, we get 



oo 



q=N+l T T 

The terms in the last expression are now estimated using a telescoping sum: 



b N 



2 ™ / >• \ 2 / * x 2 



E 



3r+2 



2iV + l/ ^\2r + l V2(r + 2) + l 

/ r=N \ / \ v ' , 



oo 



E, b r b r+ 2 \ I b r b r+ 2 
9<r 4- 1 ~~ 9.(r 4-914-1 9r 4- 1 



r=JV 
oo 



2r + 1 2(r + 2) + 1 M 2r + 1 2(r + 2) + 1 

b r b r+ 2 \ 



< 




2(r + l) + l 



< IMU^-l^lKlU^-l,!)! 



where we have used ||u||| 2 (-7-n = ]C^Lo l^</l 2 2^FI anc ^ IKIIiaer 1 ) = S^olM 2 ^!' We therefore 
conclude |(u — Hnu)(1)\ 2 < [|u[|j^a ||it'||^2. In particular, the above developments provide a simple 
proof of Lemma 13.61 below for the special case a = 0. This lemma is at the heart of the multiplicative 
structure of of Theorem 11.11 



2 Numerical results 

In this section, we illustrate the sharpness of Theorem 11.11 and Corollary 11.31 for the ID and the 2D 
case. We present the best constants in the following ID and 2D situations: 

|(ILvu)(l)| 2 < C^ tt 1 1 M 1 1 £2 (j) 1 1 M 1 1 #-!(/), 
||(IIiv«)|||2( r ) < C^lthW^^hWH^iT 2 )^ W^nuWhi^) < C^yfN + l||n|| H i (7 -2) Vu G P 2 Af(2.2) 

where I = (—1, 1) and T = (—1, 1) x {-1} c dT 2 ■ The best constants C^ lt , C^ lt are solutions of 
constrained maximization problem. For example, 

Cwt = max{ 1 1 rEiv-^ 1 1 ^2 (r) | IMI^^) |M|^i (r2) = 1, u G T 2 n}, 
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which can be solved using the technique of Lagrange multipliers. The constant C|n is more readily 
accessible as the solution of an eigenvalue problem since 



C 



2D 
H 1 



sup 

u£T 2N \\ U 



L2(r) 



Hi(T 2 ) 



The result of the ID situation are presented in Table Q] whereas the outcome of the 2D calculations 
are shown in Table [2j The 2D calculations are in agreement with the results of Theorem 11.11 and 
Cor pilar v 1 1 . 3 1 whereas the ID results illustrate |14l Lemma 4.1] and Thm. 2.2]. 



]Sf 








N 




\{n NU )(i)\* 


sup || || ,| |. 
"S-P 2 JV IMlL 2 (/)ll"llill(/) 


\\u\\ 2 


p ii ii ii ii 

u£T 2N \m\L 2 (I)\\ u \\Hl(I) 


II l|2 

v.£V 2N ll u lljjl(j) 


1 


1.1818 





8750 


55 


2.9680 


1.0462 


2 


1.8298 


1 


1436 


60 


2.9719 


1.0455 


3 


2.1527 


1 


1507 


65 


2.9750 


1.0448 


4 


2.3410 


1 


1353 


70 


2.9776 


1.0443 


5 


2.4594 


1 


1199 


75 


2.9798 


1.0438 


10 


2.7219 


1 


0826 


80 


2.9817 


1.0434 


15 


2.8221 


1 


0685 


85 


2.9833 


1.0431 


20 


2.8740 


1 


0611 


90 


2.9847 


1.0428 


25 


2.9051 


1 


0565 


95 


2.9859 


1.0425 


30 


2.9254 


1 


0534 


100 


2.9869 


1.0422 


35 


2.9394 


1 


0512 


105 


2.9879 


1.0420 


40 


2.9497 


1 


0495 


110 


2.9887 


1.0418 


45 


2.9574 


1 


0481 


115 


2.9895 


1.0416 


50 


2.9633 


1 


0471 


120 


2.9901 


1.0414 



Table 1: ID maximization problems 
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N 






l|njvtt||^ 1(T2) 


sup - — - — jp-^ 

uev 2N ll u llL 2 (r 2 )lr u llH 1 (T 2 ) 


sup — — "2 

u£V 2 N I U lljfl(7-2) 


SUp - — -5 - 


1 


1.8417 


1.4717 


0.63072 


2 


2.4820 


1.7051 


0.53460 


3 


2.8401 


1.7157 


0.50256 


4 


3.0694 


1.6988 


0.46983 


5 


3.2214 


1.6814 


0.45149 


6 


3.3282 


1.6683 


0.44284 


7 


3.4079 


1.6585 


0.43994 


8 


3.4701 


1.6508 


0.43835 


9 


3.5203 


1.6448 


0.43732 


10 


3.5619 


1.6398 


0.43664 


15 


3.6957 


1.6244 


0.43502 


20 


3.7681 


1.6165 


0.43421 


25 


3.8129 


1.6117 


0.43370 


30 


3.8430 


1.6084 


0.43346 


35 


3.8643 


1.6061 


0.43344 


40 


3.8801 


1.6043 


0.43361 


45 


3.8921 


1.6029 


0.43396 


50 


3.9015 


1.6018 


0.43445 


55 


3.9090 


1.6009 


0.43506 



Table 2: 2D maximization problems 
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3 One-dimensional results 



As mentioned in Section 11,11 many aspects of the proof of Theorem and Corollary 11.31 can be 
reduced to one-dimensional settings. In this section, we collect the univariate results for the proof of 
Theorem 11.11 in Section 13.21 and for the the alternative proof of Corollary 11.31 in Section 13.31 

3.1 Preliminaries 

We denote by Pt'^ , a, /3 > -1, n € N, the Jacobi polynomials, Prom [HI (4.3.3)] we have the 
following orthogonality relation for Jacobi polynomials and p, g 6 No: 

(1 - x) a (l - xfP^\x)P^\x)dx = ji a '^S m ; (3.1) 



here, 5 Ptq represent the Kronecker symbol and 

<a,p). l a+P+1 T{p + a + l)T{p + P + l) 

7p ' 2p + a + /3 + l p\T{p + a + P + l) ' { '> 

Furthermore, we abbreviate factors that will appear naturally in our computations: 

2(g + l) 2q + 2a 

hi{q,a) := -— — ■ — — — ■ — , gx{q,a) ■ 



{2q + a + l){2q + a + 2y (2q + a - l)(2q + a) ' 

** M := (2q + a + 2)(2q + a y 9 ^ ■= (2q + a - 2)(2q + a) ' ^ 

h ( \ 2(g + a) 2q - 2 

h 3 (q,a) := — — ■ — — - — ■ — -, g 3 {q, a) ~ 



{2q + a + \){2q + a)' ydW ' (2q + a - l)(2q + a - 2) ' 

By direct calculation we can establish relations between hi and gi. 

Lemma 3.1. Let h\, ha, h 3 and g\, g 2 , g 3 be defined in 113.3}) . Then there holds for any q > 1 and 

a £ N 



, (3.4) 



gi(q+ I, a) _ h 3 (q + l,a) ff 2 (g+ _ h 2 (q,a) g 3 (q+ I, a) _ hi{q - l,a) 

(a,0) ~ (a,0) ' (a,0) ~ (a,0) ' (a,0) ~ (a,0) 

iq ig+i Ig lq Ig lq-i 

("^-f^o) M?, «) + (- 1 ) 9+1 ^oT /l 2(g + 1, a) + (-l)^ 2 -^/^ + 2, a) = 0. (3.5) 

Tq ' Tq+1 %+2 

Furthermore, for any q > 

h 2 (q,a) - hi(q,a) = h 3 (q,a). (3.6) 

Proof. This follows directly by simple calculation and the definition of the terms. Details can be found 
in Appendix [Bj □ 

We will denote by Pq a ' ^ the antiderivative of Pq W f , i.e., 

P^°\x): = J X ^ Q _?\t)dt. (3.7) 

The following lemma states important relations between Jacobi polynomials, their derivatives, and 
their antiderivatives. 

Lemma 3.2. Let a 6 No and hi, gi, i 6 {1, 2, 3}, be given by 113. 3\) and by 13. 2i) . Then we have 



(i) for g > 1 

X (l-t) a Pf>°\t)dt = -{l-xr(h 1 (q,a)P$;\x)+h 2 (q,a)pWH 
(it) for q > 2 



P~t 0) (x) = gi(q,a)P^°\x) + g 2 (q,a)p£f\x) +g 3 (q,a) p£' 2 0) (z), 



(mj /or g > 1 



7g 7 g -i 7<? 7 g +i 

Proof The proof of (ji]) relies on relations satisfied by Jacobi polynomials; see Appendix [B] for details, 
(jn]) is taken from [2]; (pH|) is obtained by differentiating (Jn]) and using Lemma 13.11 □ 

3.2 The relation between expansions of a function and its derivative 

The essential ingredient of the one-dimensional analysis in Thm. 2.2], [171 Lemma 3.5], [141 
Lemma 4.1] and as illustrated in |(IV)| of Section 11.11 is the ability to relate the expansion coeffi- 
cients (u q ) q ^ of the Legendre expansion u = J2 q u q P q °'°^ to the expansion coefficients (5 9 )£2_ °f the 

Legendre expansion v! = ^2 q b q P q °'°^ . This relation generalizes to the case of expansions in Jacobi 
polynomials. A first result in this direction is (see also [6} Lemma 2.1] and Lemma 2.2]): 

Lemma 3.3. Let a G N . Let U G C\-l,l) and let (1 - x) a U{x) as well as (1 - x) a+1 U'(x) be 
integrable. Furthermore, assume lim(l - x) 1+a U(x) = and lim (1 + x)U(x) = 0. Then the 

x-*l x— )•— 1 

expansion coefficients 

u q := J 1 (1 - x) a U{x)p( a ' 0) (x)dx, 
b q := C (1 - x) a U'(x)Pjj a ' 0) {x)dx. 



satisfy the following connection formula for q > 1: 

u q = hi(q,a)b q+ i + h 2 (q,a)b q + h 3 (q, a)b q -i. 

Proof. Follows from an integration by parts and the representation of antiderivatives of Jacobi polyno- 
mials in terms of Jacobi polynomials given in Lemma 13.21 (i). We refer to Appendix|B]for details. □ 

This connection formula between the coefficients u q and b q allows us to bound a weighted sum of 
the coefficients u q by a weighted sum of the coefficients b q : 

Lemma 3.4. Let U G C 1 (— 1, 1) and assume j\ \U{x)\ 2 {l — x) a dx < oo as well as j\ \U'{x)\ 2 (l — 
x) a dx < oo. Let u q and b q be defined as in Lemma 13.31 Then there exist constants C\, C 2 > 
independent of a and U such that 

oo _. oo ^ „i 

g=l 79 g =0 Tq J - 1 

Proof The result follows from the relation between u q and b q given in Lemma 13.31 and from bounds 
for h\, h 2 , h$. □ 
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The following simple lemma is merely needed for the proof of Lemma 13.61 below. 



Lemma 3.5. Let a E No and q > 1. Then there exists a constant C > independent of q and a such 
that 

N 

a V -77 < 2— — \/N>q + a. 

Proof. The proof follows by the standard argument of majorizing the sum by an integral. □ 

While Lemma 13.31 shows that the coefficients u q can be expressed as a short linear combination of the 
coefficients b q (a maximum of 3 coefficients suffices), the converse is not so easy. The following lemma 
may be regarded as a weak converse of Lemma 13.31 since it allows us to bound the coefficients b q in 
terms of the coefficients u q and weighted sums of the coefficients b q . This is the main result of this 
section and the key ingredient of the proof of Theorem 11.11 as it is responsible for the multiplicative 
structure of the bound in Theorem 11.11 In the proof of Lemma 13. 6( the reader will recognize several 
arguments from |(IV)| of Section 11.11 

Lemma 3.6. Assume the hypotheses of Lemma \ 3.3l Let q £ No. Let u q and b q be defined as in 
Lemma 13.31 Then for q > 1 there exists a constant C > independent of q and a such that 



I I \ 1/2 ( \ 1/2 

V><? "i / \3>?-i 'j 



\b q -i\ 2 + \b q \ 2 < C2 a+1 



Proof. We may assume that the right-hand side of the estimate in the lemma is finite. In view of the 
sign properties of hi, hi, /13 and (|3.6p we have 

\hi(q,a)\ + \h 2 (q,a)\ = \h 3 (q,a)\. (3.8) 

We introduce the abbreviations 

h 2 {q,a) a(2q + a + 1) 



h 3 (q,a) (2q + a + 2)(q + a)' 

a(q + 2)(2q + a + 1) 

a q\ l ~ a q+l) ~ 



q " ^ ~ q ^' (2q + a + 4)(g + 1 + a)(q + a) 
By rearranging terms in Lemma 13.31 and using the triangle inequality we get 

\h 3 (q,a)\ \b g -i\ < \u q \ + \h 2 (q,a)\ \b q \ + |^i(g,a)| |6 g+ i| 

We set 



' \h 3 (q, a) I 

and by applying ()3.8p we arrive at 

< z q + a q \b q \ + (1 - a q )\b q+1 \. (3.9) 

Iterating ()3.9[) once gives 

\bg-l\ <Z q + OLg{Zq + \ + a q+1 \b q+ i\ + (1 - a q+1 ) \ b q+2 \ ) + (1 - a q )\b q+ i\ 
< z q + a q z q+ i + (l - a q (l - a q+ i))\b q+ i\ + a q (l - a q+ i)\b q+2 \ 

= Z q + a q Z q+ i + (1 - e q )\b q+ i\ + £q\bq +2 \. 
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Squaring and Cauchy-Schwarz yields 

< (z q + a q z q+ i) 2 + 2{z q + a q z q+1 )((l - e q )\b q+1 \ + e q \b q+2 \) 

+ (1 - e q ) 2 b 2 q+1 + e 2 q b 2 q+2 + 2e q (l - e q )\b q+1 \ \b q+2 \ 

< {z q + a q z q+1 ) 2 + 2{z q + a q z q+1 )((l - e,)|6,+i| + e q \b q+2 \) 

+ ((1 - e q ) 2 + e q (l - e q ))b 2 q+l + (e 2 q + e q (l - e q ))b 2 q+2 . 

If we abbreviate for the first two addends 

f q ■= {z q + a q z q+1 ) 2 + 2(z q + a q z q+1 )((l - e q )\b q+1 \ + e q \b q+2 \) (3.10) 

we obtain 

b 2 q ^<fq + (l-eq)b 2 q+l + e q b 2 q+2 , 

which we rewrite as 

b 2 q -i ~ b 2 q+1 < f q + e q (b 2 q+2 - b 2 q+1 ) . 

Next, we want to employ a telescoping sum. Since we assume that the sums on the right-hand side of 
the statement of this lemma are finite, i.e., 

E^of4 <oo < E^o) 6 ' <o °' ( 3 - n ) 

3 7? 3 7? 

and since ^ 0) < (j + a)2~ a we have y/q \b q \ — > for q — > oo. Hence, we can write 

^ 3 

OO 

+ b 2 q = ^2 b 2 q _ 1+2j - b 2 _ 1+2j+2 + b q+2j - b 2 q+2j+2 

3=0 

oo 

+2j + £ q+2j (b 2 q+2+2 j ~ bg +1+2 j) + fq+l+2j + £ q+l+2j (&g+3+2j ~ b l+2+2j) 



j=0 






oo 


oo oo 


oo 


= E fl+3 


~ E £ <?+2i fo g+l+2i + E ( £, ?+ 2 i ~ £ 9+2j+l) 


b 2 +2+2j + E g q+l+2j kq+3+27 


3=0 


j=0 j=0 


3=0 


oo 


oo oo 


oo 


= E fi+i 


- £qb 2 q+ i ~ ^ £ 9+2+2i^+3+2i + E ( £ <l+ 2 3 


- £ q+2j+l) b 2 q+2+2j + ^2 e 9+l+2j^ + 3+2j 


3=0 


3=0 3=0 


j=0 


oo 


oo 


oo 


= E fi+i 


- £qb 2 q+ i + ^ ( £ q+l+2j ~ £q+2+2j) & g+3+2j 


+ E ( £ 1+ 2 3 ~ £ 9+2j+l) b 2 q+2+2 j 


3=0 


3=0 


3=0 


oo 

= E fl+3 


oo 

— £qb 2 q+ i + ^2 ( £ Q+3 ~ £ q+j+i) b 2 q + 2 +j- 




3=0 


3=0 





We conclude, noting that e q > 0, 

b\_ x + b 2 q < b 2 q _ x + b 2 q + e q b 2 q+1 <F q + S q+2 , (3.12) 

where 

F q :=^2fj, (3-13) 

j>q 

S q := ^e'jb 2 with e) := \ £j - 2 - £j-i|. (3.14) 

i>9 
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By positivity of e'- and fj we have S q +i < S q as well as F q+ \ < F„. Therefore, we get from (|3.12p and 



3 J 3 " ovc ^g+i — 0,0 vvcli 0,0 x <?+i — ± q 
the definition of S q 

S q = £ q b q + £^+1^+1 + £^+2 

< 5 9+2 + max{e^,e / g+1 }S , g +3 + max{e' 9 , e' ?+1 }F ?+ i 

< (1 + max{E q ,e' q+1 })S q+ 2 + max{e^, 

Applying the notation 

e' q := max{e' g , e^ +1 } 

we have 

S 9 < (l + <)S g+2 + <% (3.15) 

Iterating ()3. 15|) iV times leads to 

TV TV j-1 

Sq < S q+ 2N+2 \\{^ 

j=0 j=0 i=0 



A calculation shows 



«5 £ T~T"4f = (3-17) 



From the definition of S q in (|3,14|) . (|3.11|) . and (|3.17|) it follows that lim^oo S q = 0. Furthermore, we 
can bound the product uniformly in N : 

N I N \ / N \ 

+ <+2i) = exp [J2 + < +2i ) < exp J] 4 +2j , (3.18) 

j=0 \j=0 J \j=0 J 

where in the last estimate we used the fact that ln(l + x) < x for x > 0. From (|3.17p we get 

N N N 

where we have used Lemma 1331 in the last step. Since < 1, inserting (|3.19p in (|3.18p gives 

N 

n( 1 + <+2,)<C (3.20) 

3=0 

Now, by passing to the limit N — > oo in (|3.16p . we obtain a closed form bound for S q : 

oo j—1 

S q < X) 4+2^^+23 II^ 1 + £ 9+2i)- 
j=0 i=0 

Applying (|3.19|) . (|3.2Up . (|3.17p . and the definition of F q we can simplify 

x i 

S<1 ~ Z> £ q+2jFq+2j ~ / , / , fil 1 — ^2 = / , fi / , 7 1 — ^? ~ i F q- 

{pi +3) -^{a+3) a + q 



3=0 j>q i>j i>q 3=q 
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Inserting this estimate in (|3. 12|) and using Q+ ° +2 < 1, we arrive at 

b 2 q -i +b 2 q <F q + —?—F q+2 <F q + F q+2 < F q . (3.21) 
ft T y t ^ 

We are left with estimating F q . By the definition of F q in ()3.13j) and the definition of f q in ()3.10f) we 
have 

F i = E^ + «^i+i) 2 + 2 E^ + a i^+0((l " £j)\bj+i\ + ej\bj+2\). (3.22) 

j>q j>q 

Now we estimate both sums separately starting with the first one: 

E {zj + a i z i+if % E z 1 + a l £ E 4 ( 3 - 23 ) 

i>9 i><2 j>Q 



Next, we use the relation between u q and b q from Lemma 13.31 Furthermore, we note that hs(q,a) 



> 



2 ' ; . Hence, we obtain 

2 



2 _ \ u g\ < ptt+i rfgl \ u a 

1 ~ \h„fn rvM2 ~ - 



|fc 3 (<7, «)| 2 ~ 7 («.0)fc3( g>a ) 
|u I 1 

= 2 " +1 (a!o) ^ 3 ( g a ) l^i(?' a ) b 9+i + h 2 (q,a)b q + 7t 3 (g, a)6 g _i| 
< 2" +1 ^U(l - a,)|6 9+1 | + a,|6 g | + \b q ^\). 

Tq ' 

Inserting this in the bound (|3.23p . we get by applying the Cauchy-Schwarz inequality for sums 

1/2 / \ 1/2 



We continue by estimating the second sum in (|3.22p . Using again z q < 2 a+1 \u q \/"f q a '°' > we get 
Y^izj + a jZj+1 )({l - ej)\b j+1 \ + Ej\b j+2 \) 



j>q 



^ 2a+1 E^o)(N+ a i K+il)(( 1 -^)l^+il+ ej \b j+ 2\) 

S 2 " +1 E ^0) (N + l«i+il) (to+il + to+al) 

^ +1 (e4o>w 2 ) (e^otW 2 ) ■ 

\i>,Tj ) \i>,-irj ) 

In view of (|3,2ip the last two estimates conclude the proof. □ 
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3.3 Stability of truncated Jacobi expansions 

In this section, we study the stability of the operator that effects the truncation of an expansion in 
Jacobi polynomials. We analyze this operator in the weighted P^-norm. In other word: The main 
result of this section, Lemma 13.81 generalizes [5j Thm. 2.2], where the case a. = is studied, which 
corresponds to the analysis of the ^stability of the L 2 -projection. This section is closely tied to an 
alternative proof of Corollary 11.31 and not immediately required for the proof of Theorem 11.11 

Lemma 3.7. For all a, q G No 

J (1 - x) a | (P g (a ' 0) )'(x)| 2 dx < Aq{q + 1 + af lq a ' 0) . (3.24) 

Proof. We abbreviate P q := P q a ' ^ and I q := /^(l — x) a \P' q {x)\ 2 dx. The assertion in case of q = is 
trivial. For a = see [TJ (5.3)]. A direct calculation shows 

j2 = Q 2= (« + 2) 2 2*+ 1 2= (3 + q)(a + 2) +1- 

1 4 a + V 2 2(a + l) 

the assertion of the lemma is therefore true for q G {0, 1, 2} and all a. Thus, we may assume a > 1, 
q > 2. From Lemma 13.21 (ii) with g + 1 and g there we get 



9i(g + l, a) 5i(g + l,a)ffi(g,a) 



where 

# 2 (g + 1, 0)53(9, a) a(2g + l + a)(g- 1) 



-<? •- 



51 (g + 1, a)s>i(g, a) (g + 1 + a)(2g + a - 2)(g + a) 



We note that < e q < 1. Furthermore, we calculate 



((1 - e^P^ - e 9 P^ 2 ) 2 = (1 - e q f{P' q ^f + e 2 q {P' q _ 2 f - 2e q (l - e q )P q ^P q ^ 2 
so that by integration, Cauchy-Schwarz, and < e q < 1: 

/_V - xT ((1 " " e q P' q _ 2 {x)f dx < ((1 - + e 9 V 2 ) 2 . 

By the orthogonality properties of the Jacobi polynomials, we conclude in view of (|3.25|) 



/2 = ( 1 ^ >,0) , f g 2 (g + l,«) \ (a,0) 

+ ^ (1 - *Y ((1 - ^)^-i(^) - e q P' q _ 2 (x)f dx 



^ ( 1 ^ (a,0) , 52(9 + 1, a) \ ( a ) , - r , r x2 

\g 1 {q+l,a)J q \9i{q+l,a)gi(q,a) J q 

We proceed now by an induction argument on q for fixed a. The induction hypothesis and the fact 
that q 1 y q(q + 1 + a)^' ^ is monotone increasing in q provides 

((1 - e q )I q . x + e q I q - 2 ? < 4(g - l)((g - 1) + a + 1) 2 7 ^' 1 0) , 
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and we obtain by some tedious estimates for the other two terms appearing in the bound of I 2 +1 (see 
Appendix IB]) : 



51(9 + l,a) 



,* ( * + W J 7fji 0) +%-l)(, + a)^ 
5i(9 + l,a)5i(9,a) ' ' 



^^g + lKCg + lJ + l + a)^? 



I 4 (5 + «) 2 (5-l) 7 n 



4((g + l) + l + a)V)& 0) 



4(g + l) ' 4(g + l) ' {q + 2 + Q )2 (? + i )7 £.j» 
(9 + a) 2 (9- l)(2g + a + 3) 



+ 



4((9 + l) + l + a) 2 97S'i 0) 



4(9 + 1) (9 + 2 + a) 2 (9 + l)(29 + a-l) 

2 



+ 1 



9 + « + 2 



1 



9 + 1 



1 



4(9 + 1) 

The proof is completed by observing that the expression in brackets is bounded by 1. 
We now study the stability of truncating a Jacobi expansion. 



29 + a - 1 



□ 



Lemma 3.8. Let a E No- Let u q and b q be defined as in Lemma \3.<A Then there exists a constant 
C > (which is explicitly available from the proof) independent of a and N such that for every N E N 
we have 

2 

/ a-*)" £i^(^ 0) )'(*) ^^^E tin 2 

g=0 lq q=0Tq 



Proof. We abbreviate P„ := P^ a ' % We compute 



I 1 

Y ~(^fi) U i P 'i = Y ~Koy [hi(q,a)b q+ i + h 2 (q,a)b q + h 3 (q,a)bg- 1 ]Pg 



q>N+l lq 



q>N+l lq 

Y 6 <?^7jy /l 3(9 + l,a) J Pg+i + Y b i~^0) h2 ^ a ^ P c 

q>N lq+1 q>N+l lq 

+ Y h q-^fi) h ^1- l i a ) p 'q-i 

q>N+2 lq-1 



Y b « 

q>N+2 



1 11 

-7-^^1(9 - 1, a)-P,_i + -7^5T /l 2(9, a)i^ + -7^7^3(9 + 1, a)P q+1 

lq-1 Tq ' lq+1 

N+l 



+ ~J^fi) hs ( N + 1 > a ) P N+i h N + £ ^7^2(9, a)P' q b q 
1n+i 



(a,0) 
<?=7V lq 



With Lemma 13.21 (iii) we therefore conclude 



7V+1 



Y -^ u i p q= Y b ^^ P ^^ h ^ N + 1 ' a ) P N + ibN+J2^ h ^ a K b r 

q>N+l lq q>N+2 lq 1n+1 q=N lq 
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Inserting the result of Lemma 13,71 gives 



[1-xT 



(a,o) u i P 'q 



q>N+l 1q 



q>N+2 lq 



(a,0) 

7jv 



AT+l 



CiV E^0)( iV + a ) 2/l 2(^«) 



=iv7? 



2^2 



< 



(q,0) 9" 



g>7V lq 

This allows us to conclude the argument since (Pq )'(sc) = 0. 

Lemma 3.9 (Hardy inequality). For /3 > — 1 and J7 G C^O, 1) n C((0, 1]) £/iere /joWs 



1 a^ 9 |Z7(a:)| 2 dac < (^yl / V +2 | £/» | *<te + 



2 lrr//„\l 2 j„ , 1 irr/iM2 



+ 1 



\u(i)\ z . 



□ 



Proof. This variant of the Hardy inequality can be shown using [16} Thm. 330]. See Appendix iBl for 
details. □ 

4 Properties of expansions on the tetrahedron 

To save space we will sometimes denote points in M 3 by just one letter, i.e. £ = (6,6,6) f° r points 
in 7~ 3 and 7/ = (7/1,772)%) for points in 5 3 . 



4.1 Duffy- Transformation 

We recall the definition of the reference triangle, tetrahedron, and the <i-dimensional hyper cube in 
The 3D-Duffy transformation D : 5 3 — > 7~ 3 , [TT], is given by 



with inverse 



(1 + m ){l -r? 2 )(l-r? 3 ) (l + r/ 2 )(l-r/ 3 ) 
4 1, j 1,773 



(4.1) 



•D x (6,6,6) = (viiV2,V3) 



1 + 6 1 + 6 
-2 - 1,2 - 1,6 



6 + 6 



1-6 



Lemma 4.1. TTie Duffy transformation is a bijection between the (open) cube 5 3 and the (open) 
tetrahedron 7~ 3 . Additionally, 



56 

dm 



1 3 



i(l- 77 2 )(l-r7 3 ) 











i(l + 7/0(1 -7/3) |(1 -773) 
-1(1+7/0(1 -7/ 2 ) -|(l + 77 2 ) 1 



(l- 77 2 )(l-7/3) 



4 2(1 + 7/0 2(1 
2(1-%) 1- 

(l-7/ 2 )(l-7/ 3 ) 



-m) 



det D' 



1 - V2 \ fi-m 
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Proof. See, for example, [18]. □ 

Lemma 4.2. Let D be the Duffy transformation and Y := 7~ 2 x {—1}. Then D(T) = T and D is an 
isometric isomorphism with respect to the I? (T) -norm, i.e., for sufficiently smooth functions u, we 
have for u = uo D the relation ||ii||x 2 (r) = ll^llL 2 (r) • 

Proof. Follows by inspection. □ 

4.2 Orthogonal polynomials on tetrahedra 

In terms of Jacobi polynomials Pji*'^ we introduce orthogonal polynomials on the reference tetrahe- 
dron T 3 often associated with the names of Dubiner or Koornwinder, |10lil8[ll9l: 



Lemma 4.3 (orthogonal polynomials on T 3 ). For p, q, r £ No set ipp,q,r '■= ^p,<?,r ° D 1 , where 4>p >q> ', 
is defined by 



Then the functions ip p , qi r are L 2 (T 3 ) orthogonal, satisfy i^ Pl q, r £ 'Pp+q+r{7~ 3 ) , and 

2 2 2 

VV,g,r(£)Vy,g',r'(£) d£ = 5 Pj p>5 q!q /8 r 



^ TyM „ w ry ,„ „ w - HM r > r 2p + 1 2p + 2q + 2 2p + 2q + 2r + 3 

(2p+l,0) (2p+2 9 +2,0) 
"p,p'V q , q >"r,r' Jp 22p+1 2 2p+2g+2 ' 

Proof. The proof can be found in Appendix [Bj □ 

4.3 Expansion in terms of i\) v ^ T 

Since the functions (i^ p ,q,r)p,q,re^o f° rm a se t of orthogonal polynomials, any it 6 L 2 (T 3 ) can be 
expanded as 

_ ~ (iP p , q , r ,u) LHT3) t _ ~ i 2^ 2 2 ?+ 2 ^ 2 

p,q,r=0 II r P>y' f IIL 2 (7' 4 ) p,<?,r=0 7p 7? 7r 

Wp l(? ,r := (lpp,q,r,u)L 2 (T 3 )- ( 4 - 3 ) 

A basic ingredient of the proofs of Theorem 11.11 and Corollary 11.31 is the reduction of the analysis to 
one-dimensional settings, for which we have provided the necessary results in Section As already 



flagged in (III) of Section 11.11 the ^-variable plays a special role. This is captured in Definition [5] 
below, where the functions U p>q and U p q (with expansion coefficients u p ,g, r and u' p q r ) are introduced. 
Before that, we introduce for a function u defined in T 3 the transformed function u := u o D and get 

r3 u(0^ P ,q,r(0 d£ = J u{ri)ipp tqtr (v) ( ~ 2 V2 ^j f ~ 2 V3 ^j dr] 

C / 1 \ P+l / -i \ P+9+2 
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Definition 4.4. Let p, q E No and u € L 2 (T" 3 ). We define the functions U p>q : R — > R and C/p g : 



as weZZ as i/ie coefficients u p ^ r and u' p q r by 

u P , q (m) ■= 1^ u(v)Pi m (vi)P^ p+1 '°\m) O^Y* (4.5) 
1 (i - m) 2p+2q+ %Am)P? p+2q+2 '°Hm) dm, (4.7) 



p,q,r ■— / V- 1 - 'HI ^V,1\ 



Pi?!'* 



I' (i - m) 2p+2q+2 u p Jm)P^ 2p+2q+m (m) dm- (4.8) 

With this notation, we have by comparing (|4.4p with (|4.7p 



«M,r = ^7+2 £(! " %) 2p+29+2 ^, g (r?3)P, {2p+2<?+2 ' 0) (%) d V3 = (4.9) 



Since for sufficiently smooth functions u the transformed function u is constant on m = 1> orthogo- 
nalities of the Jacobi polynomials give us 

U p>q (l) = for (p,g)^(0,0). (4.10) 
4.4 Properties of the univariate functions C/ M and [/ M 

We start with some preliminary considerations regarding estimates for partial derivatives of the trans- 
formed function u. We have 

d m u( V ) = {1 ~ m) ^~ r]3) (d 1 u) o D^), (4.11) 



^ ran) - (1 + 7?l) 4 (1 m) (d lU ) o z>fa) + ^-^(d 2 „ i o d ( a 1 2) 



Vfa) = - (1 + ??l) 4 (1 — (9m) o Ofa) - ^±^W) o Dfa) + (d 3 u) o D( V ), (4.13) 
where di is the partial derivative with respect to the i-th argument. In particular, we get 

L K%)|2 (t^) dr] = L l(diu) ° D{r])? (~?) (^) 2 dv 

= \\d^u\\h(r^ < l|Vn||| 2{r3) , (4.14) 

/ 53 l^^)| 2 (l^)^< /^K^o^)! 2 ^) 2 ^) (^) 2 ^ + ll%-llW3) 

;$ Il%«ll! 2 (r3) + H%«llL»(r>) ~ H Vu lliw ( 4 - 15 ) 

These estimates are useful to prove the following lemmas. 

Lemma 4.5 (properties of U p ^ q ). Let u G H l {T) and U p ^ q be defined in Definition \4-4\ Then there 
exists a constant C > independent of u such that 

00 1 1 ( ^ / 1 \ 2 

E (o,o) ( 2 p+i,o) / , \ U pM\ 2 I — ^ ) *?3 = llnll^^), (4.16) 
p, g =0 7p 7g ■ y_1 v 7 

£ ^^0) / , I^MI 2 HP <*& < C||V.||i 2(r3) , (4.17) 



m=0 lp lq 

oo 



E (oo) (2 P+ io) ^ + g) 2 / l^(%)| 2 ^3 < C||V U ||| 2(r3) . (4.18) 
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Furthermore, we have for F = T 2 x {— 1} 

00 1 1 

E (o,o) ( 2 P+ i,o) l^(- 1 )l 2 = IHlW ( 4 - 19 ) 



v- 



=0 7p 7g 



Proof. Proof of flUBD and ()47T7|> : The fact that P^°' 0) (r?i)Pg {2p+1 ' 0) (r? 2 ) ( ^) are orthogonal 



2 

polynomials in a weighted L 2 -space on S 2 and the definition of U p>q imply (for fixed 773) the represen 
tation 

'«<»> - E 4o ^^<^ ' 0, <'<'> F « ,W ' 0) M (^ 

Pt q=0 Jp 1q \ 

which in turn gives 

/ Hv)\ 2 0^) dmdm = E l^fe)l 2 • (4.20) 

Since detD' = ( 1 ^' 2 ^ f 1 7 2 ??3 ^ , multiplication with ^ ^ and integration in 773 gives (|4.16p . 
Similar to the representation of u above, we get for d m u 

9 V3 U(V) = E 
Reasoning as in the case of (|4.16|) yields 

E (0,0) (2p+i,o) / , l^,^)! 2 (— 2^) d % < l|Vn|| 2 2(r3) , 

p,5=0 lp 7q J - 1 \ ' 

which immediately leads to ()4.17p . 

Proof of fl4.18p : By definition, we have 

£W%) = ^r / n(^)Pf° ) ( ? ?i) J P ( ! 2p+1 ' 0) (r ?2 )(l-r ?2 r 1 dr ?1 dr ?2 . (4.21) 

We consider the integration in ?? 2 . Integration by parts then yields 

u(r ? )P^+ 1 '°)(7 ?2 )(l - r ?2 ) p+1 dr ?2 

>(r?)(l - ifc^P^f ^(ffc)) 1 - T ^ (fifa)(l - m f^) P^ lfi \ m )d m 



-1 



'-1 

Hence, we obtain by inserting into (|4.2ip 



d m u( v )(i - m ) p+l P^i m {m) + (p + - %) p ^5 +1 ' 0) (%)d%. 



s(r?) ^ o,o)(?7i)(i - ^) p n ( ?i +i,o) fe)^i^ 

^ / ^ii(i7)P^(m)(i - mT+'P^^^dmdm 



s 2 



|±i / - m) p P^i m {m)dmdm, 



18 



where in the last equation we used integration by parts in rji and the fact that Pp (t)dt = for 
p > 1. With the abbreviation gi := gi{q + 1, 2p + 1), i = 1, 2, 3 we have by Lemma [3.21 (ii) for p, q > 1 
the following relationships: 

Pf p+lfi) (m) = 9iP^ +m (v2) + 9 2P^ +1 '°\m) + 53n (2p i +1 ' 0) M> ( 4 - 22 ) 

3°' 0) foi) = ^7 (iftV) " • (4-23) 
Furthermore, we introduce two abbreviations 

p+i 



ZpM ■= Jj9 m u)(ri)P^ H m ) Pt +m (V2)dmdr)2 (4.24) 

ZpM ■= j s {{dm) o D){ V )P^\ m ) {^-J^ Pt +m (V2)dvid V2 . (4.25) 
Since we have KW\ . using (ET221) . (Q3|) . (ET24"j) . and EH| we get 



Up, q {m) = -{9lZp,q+l(V3) + 92Z p , q (V3) + g3Zp,q-l(m)) 



+ 



p+l\ ( 1 - 773 



0l(%H,g+l(f&) - Vl,«+l(l3)) 



2p + l 

+ 52 (ffc) - Zp-l,q(.m)) + 53 (^p+l.g-l (%) ~ VU-lt^)) 



We use pi, 52, 53 < ^ to arrive at 



2 /l - \ 

(P + < E ( 7"^ ) (\Zp+hq+l-j(m)\ + + |2p,g+l-j(%)|- (4.26) 

To estimate the terms on the right-hande side we note that the abbreviations z pq and z Pt q lead us to 
the representations 

oo 



p,q=0 lp lq ^ ' 



p, g =o 7p 7g 

Since the polynomials Pp°' \rji) ( 1 ~ 2 r?2 ~) Pg 2p+1 '°^ (772) are orthogonal polynomials on the reference 



2 

triangle T 2 we have 



/ \d m u(v)\ 2 ( V^) d ^ 2 = £ -m-wm ^ 3)|2 (427) 

K^fi) o D)(r,)| 2 *7icfc» = E (o,o) ft^ WI 2 - ( 4 - 28 ) 



p , g =0 7p 7g 

Formula (|4.27p together with an integration in 7/3 and an application of (I4.15P gives 

E (0,0) ( 2p +i,o) / \ z vM\ 2 dm = f \d m u(v)\ 2 0-^) dr > ~ ll Vu lli 2 (r«)- 
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From the representation (|4.28|> we get by a multiplication with 
use of QU1J) that 



2 



, an integration in 773 and the 



E 



1 1 



(0,0) (2p+l,0) 
p,q =0 7p 7i 



1 

-1 



z P ,q{m)\ 2 



1 - m 



dm 



\((d lU ) o D)( V )f 



1 - m \ (i-m 



dv < \\Vuf L2 



(T 3 )- 



The last two results in combination with (|4.26p conclude the proof of (|4.18p . 

Proof of (I4.19|) : For the estimate (|4.19|) . we use (|4.2Up with 773 = —1. Noting Lemma HJ we 
have 



I II 2 — II ~ l|2 

l n llL 2 (r) — H n llL2( r ) 



|«(*7i>»72,-l)| c%d?72 



and therefore the result follows. 



□ 



The estimate (|4.18|) does not provide a bound for the special cases p = or q = 0. The treatment of 
these two special cases is the purpose of the following lemma. 

Lemma 4.6 (Properties of J7o,g and U Pi o). Assume the same hypotheses as in Lemma \4-5\ Then there 
exists a constant C > independent of u such that 



00 1 1 f 1 

Yl nm^m q2 / i^^)! 2 d ^ < c\\v u \\ 2 L2{T3) 

q=0 % r q J - 1 



00 1 1 f 1 

E (0,0) (2 P+ i,o) ^ 2 / t KAm)\ 2 dm < c\\Vu\\ 2 L2{T3) . 

n l\ J-l 



(4.29) 
(4.30) 



P =o Iv 7o 

Proof. The proof of (|4.29p start with the observation 



U, 



0,9 



1 / /•! 



-1 V-l 



1 -»72 



*?2 



We expand for fixed 773 the function 772 1— > f —1 u{rj\, 772, 773)^771 as well as its derivative in terms of 

Pq . Lemma l3~4l then yields 

00 „i 00 „i 

T,nm^m q l lUoAmtdm^^^^iq + i) 2 \u , q (m)\ 2 dm 

q =0 7o ' 7? ' 0=1 7<j ' J ~ l 



< 



1 /•! 



-1 J -I 



-1 



d m u(rj)dr)i 



(1-772)^772^7/3 < / |^ 2 tt(77)|' 
.As 3 



1 - m\ , / MY7 ,,2 

dv ^ l|V-u|| L2(r3) , 



where we appealed to (|4.15|) in the last estimate. 



Analogously we deal with (|4.30p . where we expand the map 771 1— > J_-yu{r}x,r]2,m 



1—772 



p+1 



dm 
□ 



again for fixed 773 in terms of P p ' ^ and conclude with Lemma 13.41 and (|4.14p . 

The estimates for the functions U Pt q imply easily corresponding bounds for the functions U Piq : 

Lemma 4.7 (properties of U p<q ). Let u 6 H l (T) and U Pj q be defined in Definition \4-4\ Then there 
exists a constant C > independent of u such that 



1 ^ 1 1 

4 2^ (0,0) (ap+1,0) 

p,q=0 lp lq 



1 00 

E 



1 1 



4 ^ ^(0,0) (ap+i.0) 

p,q=0 IP lq 



1 _ m fV+W 



u' p ,q{m) dm < c||v-u|||2 (T 3 



(4.31) 
(4.32) 
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Proof. We have U m ( V3 ) = (1 - V3)~ ip+q) U p , q (rj 3 ) and therefore U p>q (r, 3 ) = (1 - + 
(p + ?)(!- r/s)-^^ 1 )^,^^). Hence, 



1 Z" 1 



4 7-i 





















2 ^</_;( 





2 



2 



73 

|^, 9 (%)| 2 dr? 3 + (p + q) 2 J^ \U M ( m )\ 2 dr) 3 . 



2 

Using the results of Lemma 14.51 concludes the argument. □ 

These results allow us to get bounds for weighted sums of the coefficients u p> q tr and u' r given in 
Definition I 



Corollary 4.8. A ssuTTic the hypotheses of LcuiTfid ^^. 7| and let up s q } r7 ^pqr 

be given by Definition ^. 4\ 

Then there exist constants C independent of u such that 

°° 1 1 1 

Yl Hofi) (2p+l,0) (2p+2g+2,0)l*W| 2 - C W U \\h(T s )> ( 433 ) 
p,q,r=0r P Tq 7r 

°° 1 1 1 

Yl (0,0) (2p+l,0) (2p+2q+2,0)K,q,r\ 2 < C\Nu\\l 2{T3 y (4.34) 
p,q,r=0 Ip lq lr 

Proof. Since the polynomials p^ 2 ' p+2q+2 ^ are orthogonal polynomials in a weighted L 2 -space, expand- 
ing U p%q yields the representation U p , q (r, 3 ) = YITLo tfffk+m ~"p,g,r-Pr (2p+2g+2 ' 0) (m) an d therefore 

f>l 00 -I 

/ (l-.3) 2 ^ +2 |^( % )|^3 = E^ 



. (2p+2g+2,0) |U P'9' r| 
=0 7r 



The statement (|4.33|) now follows directly from (|4.3ip of Lemma 14. 71 Analogously, we deal with (|4.34|) . 
where we expand U p q and conclude with (|4.32p of Lemma 14.71 □ 

4.5 Connections between u nnr and u'„ _ 

As we have mentioned in |(IV)| of Section [TTT1 the multiplicative structure of the estimate of Theorem ll.il 
is based on a relation between the coefficients u P) q )T and u' pqr . This connection is essentially a one- 
dimensional effect and follows from Lemma I3T 



Corollary 4.9. Let u p , q , r and u' r be as in Definition ^. 4\ and h\, h>2 and h 3 be given by i3.3\) . Then 
for r > 1 and p, q > there holds 



u, 



,,,, r = hi(r, 2p + 2q + 2)u' pqr+1 + h 2 (r, 2p + 2q + 2)t£ + fc 3 (r, 2p + 2q + 2)^ 9r _ 1 . 



Proof. By density we may assume that u £ C^iT 3 ). Lemma 13.31 then implies the result. For a 
detailed version of this proof see Appendix [Bj □ 



5 Trace stability (Proof of Theorem 11.11) 



We start with establishing a representation of the transformed function u = u o D on the face T 
T 2 x {-1}. Using (fOJl and we get 



p r (2p+2 ff +2,0) ( _ 1 j = (^rj. 



(0,2p+2 9 +2) 



(1) 



l) 7 
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Therefore we have for £ € T 

= <W O " 1 (a, = (-if 

Applying the expansion of u in (j4.2j) we arrive at 



l-r7 2 \ p 



^ 0) (^)n (2P+1 ' 0) fe) 



«(??i,?/2,-l) = u(£i,£ 2 ,-l) 

00 1 o2p+l 

Z> (0,0) (2p+i,o)- r P 

p,q,r=0 IP lq 



1 mY p (2p+l,0)/ Wiy.. 22P+2g+2 (m) 

^ (%)( 1J Up.g.r (2 p + 2g+2,0)- ^--U 



7, 



In view of Lemma 14.21 as well as the expansion (|4.2|) . we obtain for the L 2 (r)-norm 



i2 11 ~ n2 V ^ 

'ill r2^ = llltll T2rr\ = 



1 1 



L2(r) - H"IIl2(d - (0,0) (2p+i,o) 

p, g =0 7p 7g 



2 2p+2 9 +2 



r=0 



'P,<2,r ( 2p +2g+2,0) 

7r 



(5.2) 



Before preceding further, we recall [(TV)] of Section [TTT| where we studied (u — U.^u)(l) in equation 
(jl.7p . There, we observed that the infinite sum X^jv+i reduced to a short sum with merely two 
terms. A very similar effect takes place here in the multi-dimensional case in that the infinite sum 
over / in (|5.2j) can be expressed finite sum: 



Lemma 5.1. Let N > 1. Then there holds 



2 2p+2 g +2 



2P+1 



(2 P +2 q +2,0j u P^ = ("I) h 2 (N, 2p+2q + 2) {2p+2q+m u' p , q>N 

r=N lr 7jV 

N 



+ (-l) r+1 h 3 (r + l,2p + 2q + 2: 



r=N-l 

Proof. We abbreviate n pq := 2p + 2q + 2. In view of Corollary 14.91 we have 



(2p+2<?+2,0) P<Q>' 
7 r +l 



2 n Pi 



--N 7 



=jv 7 



2P+1 



--N J. 

00 

£ (-i r %(r-l, 

r=N+l 



+ /i 2 (r, n pq )u' p q r + h 3 (r, n pq )u' p q r _A 



2P+i 



n 



I'D (nj,„0) 
7r— 1 



+ ^(-l)^2( 



r=N 



2P+i 



lr 



r, n 



^,9,r+ E (-l) r+ ^3(r + l,n ? 



2P+9 



r=7V-l 



7r+l 



r=7V+l 



/ii(r - l,n P9 ) h 2 (r, ^pg 



(ri pq ,0) 
/j — 1 



+ 



2P+q ^ 2 p+(? 

+ (-l) N h 2 (N,n pq )^—^u' p ^ N + (-IY +I h3(r + I,n nn )^^u' 



p<i) 

(n pq ,0) 
lr 



7 r +l 



7at r=AT-l 7 r +l 

By (|3.5p . the expression in brackets vanishes and that concludes the proof. 



w („ p „o)"W 



□ 

Since Lemma 15.11 assumes N > 1, the terms corresponding to r = in (|5.2[) are not included. We 
study this special case now in Lemma 15.21 
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Lemma 5.2. Let u S H 1 ^ 5 ) and consider the representation of the norms in 115. £)) . For r = i/iere 
exists a constant C > independent of p, q, and u such that 



1 1 

E (0,0) (2p+l,0) 
p,g= 7p Tq 



2 2p+2g+2 
U P'9'° (2p+2g+2,0) 

To 



< C \\ u \\L^(T 3 )\\ u \\m(T 3 )- 



Proof. We first note that the coefficient |ito,o,o| — VT^I IMU 2 (T 3 ) so that we may focus on the sum 

with (p,q)^ (0,0). 

Since -y (0 ' 0) - -y (2p+1 ' 0) - 22p+1 and ^J 2 ^ 2 ^ 2 ^) - 2 2 ^+ 2 ^+ 3 t 

since 7 P — 2p+1 , 7g - ana 7 — 2p +2 q +3 ' we ° 

2 



2 2p+2q+2 
U P,qfi (2p+2 g +2,0) 

7o 



< 



- E (i,0) (2p 1 +l,0) ^ + g+ 1 ) 2 |^0| 2 
p, g=0 7?> ' Tq 



To bound the sum on the right-hand side, we note that an integration by parts and (|4.10p give (for 
(P,3)^(0,0)) 



(5.3) 
(5.4) 



2 p+q+2 u P , q ,o = I (i-m) p+q+2 u P , q (v3)dm 



p + q 



These two equations yield two estimates for u Pj q t o- From (|5.3p . we get with the Cauchy-Schwarz 
inequality 



\2 P+q+ \,qfi 



(l- m )^ 2 U p>q { m )d m 



< 



(Y ((i - m) p+q+1 f dm) (J\i ~ m) 2 \u P , q (m)\ 2 dm 

2 2p+2q+5 ,1 /1_^ 2 



2p + 2^ + 3 y_! 

From (|5.4p . we obtain a second estimate as follows: 

2 



|C/p, 9 (7? 3 )| 2 ^3 



(5.5) 



\2 p+q+2 u P , q , \ 2 < 



(p + q + 3) 2 



2 W6 |(/ M (-1)| 2 + 



(i _ m ) p+q+3 u p , q (m)dm 



where 



(i - m) p+q+ *u' M { m )d m 



-i 



< 



(J 1 ((i - m) p+q+2 ) 2 dm) (J\i - m) 2 \u' PA {m)\ 2 dm 



2"2p+2q+7 f 1 f\ 



2(p + q) + 5j_ 1 
Inserting this in the bound before yields 



m 



K, q (V3)\ 2 - 



\2^+ 2 u p>q n\ 2 < 2 



22p+2 9 ^ 



(p + q + 3) 2 



\U p , q (-l)\ 2 + 



1 



p + q + 2 J_ x 



1 fl-V3 



K q (m)\ 2 dm ) • (5-6) 



Next, we abbreviate 



p,q 



' ll ^r) 2 \ U vM?dm, r 2 >q :- 



1 fi-m 



\u' p , q {m)\ 2 dm- 
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Hence, applying (|5.5p and ()5.6p we have in view of the elementary observation min{a 2 ,6 2 + c 2 } < 
b 2 + min{a 2 ,c 2 } < b 2 + |a| |c| (for real a, b, c): 



\2 p+q+2 u Piqfi \ 2 < min 



r 92p+2g+5 



2 2p+2q+6 



< 



4 2 2 

-|J7 M (-1)| H — t„.„-, 



p + g + 2 



p,i'p,i 



This leads us to 



K,<?,ol 2 < 



(p + q + 3) 2 



l^, 9 (-l)| 2 + 



(5.7) 



Hence, we conclude 



11 11 11 

E (0,0) (2 P+ l,0) (P + g + 1 ) 2 |^,0| 2 < E (0,0) (2 P +1,0) I^(- 1 )| 2+ E 



n (0,0) (2p+l,0) 
p, ? =0 7p 7<? 



p,q=0 IP lq 



(0,0) (2^+1,0)^^ 

p,q=0 7p 7g 



< 



Mli^r) + IMIl 2 (T 3 )II v,u IIl 2 (t 3 )> 



where, in the last inequality, we used Cauchy-Schwarz for the sum involving a pq and T pq and appealed 
to Lemma 14,51 The proof of the lemma is now completed with the aid of the multiplicative trace 



inequality |M|| 2 ( r ) < IMIl^T 3 ) IMI^MT 3 ) ( see > e -S-, P> Thm - 1-6.6]). 



□ 



We conclude this section with the proof of Theorem 11.11 and Corollary 11.21 

Proof of Theorem 11.11 and Corollary 11.21 In view of the multiplicative trace inequality 
IHIi 2 (r) ~ IMIz 2 (T 3 ) IMIff^T 3 ) ( see > e -&-> El Thm. 1.6.6]) we will only show the statement \\u — 
njyu|| 2 2(r) < ||u|| £ 2( r 3)||tt|| fr i( T 3). We abbreviate n pq := 2p + 2q + 2 and c pg := ( * 0) ^+1,0) ■ % (JE2D, 

Ip Iq 

we have to bound 



n N u\\l2 {r) = c pi 

p,q=0 



E 
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E 

p+q<N 
N N-p 

EE 



1 "i 



E (-!) r 

r=N+l— p-q 7i 



{O.JV+l-p-g} 

2 

2 n Pi 



2 n pq 



(n Pq ,0) U P>1> r 



+ E 



pq 



p+q>N+l 



E(-D r 4 



2 n pq 



=0 7 



(n pq ,0) "P>1> r 



V-i 



E (-D r 



2 n pq 
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00 00 



r=N+l-p-q 7 
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N 00 

+ E E 
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L py 



=0 7 



E E< 
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p=JV+l q=0 



00 nil 

^— \ , 2 vq 

r=0 7r 
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AT AT-p 

EE' 
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!><-! 



2 n Pi 



r=Ar+l-p-(j 7r 



AT 00 

+E E c p<? 

p=0 g=Af+l-p 



=1 7r 



Si 



00 00 



+ E £■ 

p=7V+l q=0 



pq 



^ l J ^(n M .0) Up,, * r 
r=l 7r 



2 TV 



+E E 

p=0 q=N+l—p 



-pq 



=■■ S 2 
2 n pq 



Jn pq fl) U P^ 
10 



--■s 3 



s 4 
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oo oo 



+ E E 1 

p=N+l q=0 



I" I 



76 



s 5 



Lemma 1531 immediately gives S4+S5 < ||u||£2(7-3)||ii||#i(7-3). From Lemma I5TT1 as well as the estimates 
h 2 {N+l-p-q,n pq )<^ < 1 p + 9 = 0,...,JV 

h 3 (Ar+l-p- 9j n p? ), h 3 (jy + 2-p- gl n w )< ~^, + g < -i, p + <? = 0, . . . , iV 



and 



2N + 5 



(n p „0) 2 2 P+ 2 9+ 3 ' 
77V+l-p-<j 



2N + 7 



(n P9 ,0) 2 2 P+ 2 1+ 3 
^N+2-p-q 



we obtain for S\ 



N N-p 



SisEE 



o-(p+g+3)-/ 



p=0 5=0 

Analogously, we get for S2 and 5 3 

AT 00 



+ 



9-(P+9+3)f/ ^ 

^ w p,q,N—p—q 



^<E E | 2 ^ (p+9+3) ^,i| + | 2 ~ (p+9+3) < 9 ,o| ) 

p=0g=7V+l-p ^ ' 
00 00 / 2 2\ 

s 3 < E E^(| 2_(p+9+3) ^i| +|2- (p+9+3) < 9 ,o| • 

Applying Lemma f3.6l the powers of two in the estimates above and in the lemma annihilate each other. 
Hence, S\ + Si + S3 gives 



N N-p 

si + s 2 + s 3 <EE 

Cpq 

p=0 q=Q 



1/2 
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(n Pg ,0) '<W 



<->N+l-p-q 7 
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<->N-p—q 7 



(n p „0) '"JWl 



N oo 



1/2 



+ E E c ™ Et! 



p=0 g=JV+l— p 



r>i 7; 



V 1 , 

2^ (n p „0) 



1/2 



->o 7'. 



00 00 



1/2 



+ E E c ^ E^; 



< 



p=Ar+l q=0 

Vu\\ 



->1 7 



(n p9 ,0) 



(n Dl 



1/2 



->o 7 



u' I 2 

(n P9 ,0) i u P,9,r\ 



I u IIl 2 (t 3 )II Vw IIl 2 (t 3 )' 



where in the last estimate, we have used the Cauchy-Schwarz inequality for sums and Corollary 14.81 

Since ||Vtt| 1^2(7-3) < |M|/fi(r 3 ) this concludes the proof of (|1.2p . 

The estimate (|1.3p follows directly from ()1.2p in view of [27} Lemma 25.3]. 

Finally, f|l .4[) is obtained in a fairly routine way from the stability of just shown. Specifically, 
For arbitrary v G Pat we have by the projection property of ILv as well as the continuity of the trace 
operator 7o : ^^(T 3 ) -»• L 2 (3T 3 ) (cf., e.g., [23 Sec. 32]) 

||u - n iV M||i2(9 T 3 ) < ||u - ^11^,2(97-3) + \\U N (u - v)\\ L 2 {ari) < C\\u - v || B i/2 (r3) - 
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Hence, \\u — Hnu\\ L 2(q-j-3\ < C inf \\u — v\\ i/2,_. 3V Fix s > 1/2. Let ILv : L 2 (7~ 3 ) — > Vn be an 
approximation operator with simultaneous approximation properties in a scale of Sobolev spaces, viz., 

||u — njvii|| L 2(7-3) < CN~ s \\u\\jj s r'j-3\, \\u — IIjv"u||^s(7-3) < C||u||^f S (7-3) Vn G H S (T 3 )- 

(This can be achieved, for example, by combining the approximation results of [211 Appendix A] 
for hyper cubes with the well-known extension operator of Stein, [251 Chap. VII].) By interpo- 
lation theory, we then have that Id— 11^ is a bounded linear operator H S (T 3 ) — > b\^(J~ 3 ) = 
(L 2 (r 3 ),^(r 3 )) (s _ 1/2)/Sil with norm || Id -5^1^/^^^ < CJ\r«(-i/2)/«. □ 



6 /^-stability (Proof of Corollary [US) 



Our procedure to study the i? 1 -stability of the L 2 -projection II/v is to consider on the hyper cube 5 3 
the derivative d m VN of the transformed function vn = v^oD, where vn '■= TInu- This provides infor- 
mation about a directional derivative of on 7~ 3 . Through affine transformations of the tetrahedron 
7~ 3 , information about the full gradient of v at can be inferred. 

The key step is therefore to control d^Hiyu, where we denote H^u := (IIvu) o D. This is the 
purpose of the ensuing lemma. 

Lemma 6.1. There exists a constant C > independent of N such that 

(l-r /2 )(l-r ?3 ) 2 d m Il N u( V ) 2 d V <CN\\Vu\\ 2 L2(r3) \fu£H l (T 3 ). 

S3 ^ ; 

Proof. We abbreviate n pq := 2p + 2q + 2. We see that 
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by recalling the relation between u v ^ r and u v ^ r . Differentiating with respect to 7/3 shows us that we 
have to estimate the two terms 
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EE E (0,0) (2p+l,0) (npgfl) 
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3.2) 



First, we consider ()6.ip . From Lemma 13.81 with a = n pq we get 
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=0 It 
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where in the last step, we appealed to Corollary 14.81 Thus, we arrive at the desired bound for I\. 
Next, we consider (|6.2|) . We have 
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Lemma [3. 91 with f3 = 2p + 2q and the normalization convention for Jacobi polynomials (|A.4p now yield 
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(n pq ,0) 



U 



P,Q,r 



where we used the bound for I\ proved above and the multiplicative trace inequality [31 Thm. 1.6.6]. 

□ 

Proof of Corollary 11.3b For a function v and the transformed function v = v o D, the formula 
(|4.13|) provides a relation between d m v and Vu. Rearranging terms yields 



1 + 6 

Therefore, when transforming to 7~ 3 in Lemma 16. II we get 



d 2 v(0 + d 3 v(0- 



i-6 



< iv||Vu||| 2 (7-3 



(6.3) 



By the symmetry properties of 7~ 3 , we see that also the following two other permutations of indices 
are valid estimates: 
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l±^d 2 n NU (0 - i±ii d 3 u N u(0 + diUNutt) 
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(6.5) 



We abbreviate a(x,y) := — := a(6>6) anc ^ 
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-4(6,6,6):=! sym 1 + a| 3 + a| x a 23 + a 2 ia 3 i + a 32 

sym sym 1 + a 31 + a| 2 
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Hence, we see that by adding (|6.3|) . (|6.4p . and (|6.5[) we arrive at 

J^(VU N u) T A(OVU N u d£ < N\\Vu\\ 2 L2{r3) 
Next, we observe that near the top vertex (—1,-1,1), we have 



1 + 6 



1-6 



< 1 and 



1 + 6 



i-6 



< i. 



This implies that the functions a\3 and 023 are uniformly bounded on 7~ 3 . Analogously, we get bounds 
for 012, 032 and 021, 031 by studying the vertices (1, —1, —1) and (—1, 1, —1). Therefore, we have 

SUp ||A(£)||L°c (r 3) < 00. 

£eT 3 

By construction, the matrix is (pointwise) symmetric positive semidefinite. Our goal is to show 
that ^4(^) is in fact positive definite on the set that stays away from the face F opposite the vertex 
(—1,-1,-1). This can be done with techniques as in |13j by establishing lower bounds for the 
eigenvalues of A{£). A direct calculation reveals 

detA(e) = 16 ji + ggi6 + 266 + ggi + e 2 + 1 + 26 + 266 + 26 + 6 2 
<U " " (-i + 6) 2 (-i + 6) 2 (-i + 6) 2 

(6 + 6 + 6) 2 + 2(6 + 6 + 6) + i 1C (1 + 6 + 6 + 6) 2 



16 ^ ^ ^ — = 16 



v-i + 6) 2 (-1 + 6) 2 (-1 + 6) 2 " (-1 + 6) 2 (-1 + 6) 2 (-1 + 6) 2 ' 

The face opposite the vertex (—1,-1,-1) contains the vertices (—1 — 1,1), (—1,1,-1), (1,-1,-1) 
and is given by the equation 6+6 + 6 + 1 = 0- Furthermore, we conclude that the signed distance 
of an arbitrary point £ from this face F is given by 

dist(£,F) = ^(6+6 + 6 + 1)- 

Let, for arbitrary 5 > 0, 

T S :={££T 3 \ dist(£,F) < -5}. 
Then, since we stay away from the face F, it is clear that there exists Cg > such that 

det A(£)>C S V(eTj. 

Combining the above findings, we have that on T$ the matrix A(£) is in fact symmetric positive 
definite. Since the entries of A(£) are uniformly bounded in £, Gershgorin's circle theorem provides a 
constant C upper such that all eigenvalues of A(£) are bounded by C upper . 

A lower bound for the eigenvalues is obtained as follows: Denoting for fixed £ € T§ the eigenvalues 
< Ai < A2 < A3, we get from det A = A1A2A3 

Cg < det A = A1A2A3 < \\C\ pper . 

This provides the desired lower bound for Ax- Thus, we conclude that for every 5 > we can find 
c<5 > such that A(£) > cgl on T$. Hence, 

c s \VIL N u\ 2 dt < J^(VU N u) T A(0VU N u d£ < N\\Vu\\ 2 L2{Ta) . 

Affine transformations allow us to get analogous estimates for the sets that stay away from the other 
faces of 7~ 3 . We therefore get the desired result. □ 
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A Properties of Jacobi polynomials 

We have the following useful formulas (see [T81 P- 350 f], [25]): 
Recursion Relations 

a^SfO*) = («n + a 3 n x)Pt^(x) - aiP^i\x) (A.l) 

with 

a l n := 2(n + l)(n + a + /3 + l)(2n + a + /3) 
a 2 := (2n + a + /3 + l)(a 2 - /3 2 ) 

:= (2n + a + /3)(2n + a + /3 + l)(2n + a + /3 + 2) 
a A n := 2(n + a)(n + /3)(2n + a + /3 + 2) 



^(x)^P^(x) = b 2 n (x)P^(x) + 6^ ( X )p(«f )( X ) (A.2) 

with 

b l n {x) := (2n + a + /3)(l -x 2 ) 

6 2 (x) := n (a - /3 - (2n + a + /3)x) 

6^(x) := 2(n + a)(n + /3) 



Special Values 



P (a ' /3) = 1 (A.3) 
^(l)^^ ) (A.4) 
P(^)(-X) = (-1)"P(^)(X) (A.5) 



Special Cases For the Legendre Polynomial L n (x) there holds 

L n (x)=P<°>°\x) (A.6) 



Miscellaneous Relations 



d_ 

ch- 
in /~ (1 - t) a (l + tfP^\t) dt = -(1 - x) a+1 (l + xf +1 P^X l ' P+1 \x) (A.8) 



^-P^^(x) = \{a + /3 + n + (A.7) 



-i 
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B Details for selected proofs 



B.l Selected proofs for Section [3] 
B.l.l Extended proofs of Lemma I3TT1 

Proof of \3. 5\) . By definition of 7p"'^ we obtain in particular 

lq " 2q + q + 1 ' 
which leads in combination with the definition of hi, /12 and h 3 to 

(- 1 )"^0) M*> «) + (-l) ,+1 -7^oT M« + 1, a) + (-iy +2 ^h 3 (q + 2, a) 
Tg ' Tq+l Tg+2 

^q + a + l 2q + 2 t f +1 2g + a + 3 2a 



v 7 2 Q +! (2g + a + l)(2< ? + a + 2) v 7 2°+! (2q + a + 4)(2g + a + 2) 

+2 2g + a + 5 2(g + a + 2) 

1 j 2 Q +! (2g + a + 5)(2g + a + 4) 

^ g + 1 (2g + a + 3)a + a + 2 



2 a \2q + a + 2 (2q + a + 4)(2g + a + 2) 2g + a + 4 
_ / (g + l)(2g + a + 4) + (2g + a + 3)a - (g + a + 2)(2g + a + 2) 

2 a V {2q + a + A){2q + a + 2) 

Simply multiplying out the numerator concludes the proof regarding the first equation. 

Inserting the definition of hi, h 2 and h$ also leads in the case of the second equation to the 
conclusion 

2a(2g + a + l) + 2(< / + l)(2< 7 + a) 

h 2 {q,a) - hi(q,a) - 



{2q + a)(2q + a + l)(2q + a + 2) 

4q 2 + 4q + 6qa + 2a 2 + 4a 
(2q + a)(2q + a + l)(2q + a + 2) 

(2g + a + 2)(2q + 2a) 
(2g + a){2q + a + l)(2q + a + 2) 



h 3 (q,a). 



□ 



B.l. 2 Extended proofs of Lemma EOl 

Proof of Lemma \3.H!\ (i). Using rearranged versions of (|A.2|) . (|A.7p and ()A.8j) we obtain 

GLHJ 1 



J {i - trp^\t)dt + x)(i - x^p^c*) 



la-^a-xrp^w 



^ _(i _ x) a _L (1 _ ^ 2 d ^) (x) 

2g g + a + 1 ax y 

JX21 „ _, a l 1 q(a-(2q + a )x)Pi a '°\x)+2q(q + a)P^f\x) 



-(l-xy 



q q + a + 1 2g + a 



_ aP,P u; (s) + 2(g + a)P£'?(x) ~ gg + a)xP t 
" l Xj (g + a + l)(2g + a) 
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(|A,ip allows us now to replace the term xPq (x) by terms involving Pq"\ (x), Pq (x) and (x). 
Hence, we get 



(l-t) a p( a '°\t)dt = -{l-x) c 



1 



(q + a + l)(2q + a) 



aP^ \x) + 2(q + a)P^f ] ( 



(2q + a + l)(2q + a + 2) 



2(q + l)(q + a + l)(2q + a)p£f ) (x) 



+ 2q(q + a )(2q + a + 2)P^f\x) - (2q + a + l)a 2 P^°\x) 



Rearranging terms gives 

{l-t) a P^\t)dt = -{l-x) c 



1 \ 2(q + l)(q + a + l)(2q + a) p(a>0) 

(q + a + l)(2q + a)\ (2q + a + l)(2q + a + 2) q+l {X) 



2q + 2a + 2 ( 0) g + a + 1 (a , c 

+ a 2 9 + a + 2 ^ (Xj + 49 + aJ 2g + a + l" ^ 

-(1 " *)4 -7 r ife 0) (s) 

V 7 I (2g + a + l)(2g + a + 2) 9+1 V ; 



hi{q,a) 



+ 



2a 



(2g + a + 2)(2 9 + a; 



p(a.0) f ,^ + 2(q + a) p( a '°V^ 

^ W + (2g + a + l)(2g + a)^- 1 (Xj 



h 3 (q, a) 



□ 



B.1.3 Extended proofs of Lemma [3771 
Lemma B.l (details of the proof of Lemma 13. Tj) . Write 



q,a 



(i- x y 



P^)'(x) 



dx 



Then 



(i) fora = we have I* fi = q{q + 1) < K ■ q{q + l) 2 ^ '^ with K = 1 
(ii) for a > 1 we have 



T 2 

1 l,c 



l 2,a 



(a + 2)22 " +1 < K ■ 1 • (1 + 1 + a) 2 7 J Q ' 0) with K = 1/2 
4(a + l) 

(3 + a)(a + 2)2 a+1 ^ ^ „ ,„ , , , s 2 fa,o) 



(B.l) 



< K • 2 • (2 + 1 + a)^ * j with K = 9/16 (B.2) 



2(a + l) 

Proof. Proof of @): We have I? = g(g + 1) by [TJ (5.3)]. The statement therefore follows directly. 
Proof of $vj$: A direct calculation shows 

(a + 2) 2 2 Q+1 (3 + q) 



T 2 

1 l,a 



1 



(a + 3) 

4 a + 12 a + 1 (2 + a) 2 ~ 4(a + 1) 



l.(l + l + a )2 7 K0) 

this last function is monotone decreasing in a so that its maximum on [l,oo) is attained for a = 1, 
which has the value 1/2. Analogously, we proceed for the case q = 2. We have 

1 (a + 2)(q + 3) 2a+1 5 + a _ (a + 2)(a + 5) 



T 2 



'2- (2 + l + a) 2 7 £ 



2 o/ («.0) 



2(a + l) 



2 a + 1 2(3 + a) 2 4(a + l)(a + 3)" 
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Again, this last function is monotone decreasing on (0,oo) so that its maximum on [l,oo) is attained 
for a = 1. □ 

Lemma B.2 (details of the proof of Lemma l3.7p . Define for a, q £ No 

52(5 + 1,0)53^,0) a(2g + l + a)(g- 1) 



gi(q + l,a)gi(q,a) (q + 1 + a)(2q + a - 2)(q + a) ' 
Then, for a, q > 1 we have < e q < 1. 

Proof. Clearly, e q > 0. To see the estimate e g < 1, we have to show 

? 

a(2g + l + a)(g-l) < (g + a)(g + 1 + a)(2g + a - 2) 
a(g-l)(2g + a-2) + 3a<g-l) < (2g + a - 2)(g + a)(q + a + 1) 
<^=> (2g + a-2)(a(g-l)-(g + a + l)(g + a))+3a(g-l) < 



>i 

This last inequality is certainly true if 



(g + a 


+ l)(g + o) + 


3a(g — 1) 


? 
< 





4a(g - 


-1) - (g + a + 


l)(g + a) 


? 
< 





Aa(q - 


-l)-(g + a) 2 


- (q + a) 


? 
< 





-1)- 


q 2 — 2aq — a 2 


- (g + «) 


? 

< 





2aq - 


- 4a — q 2 — a 2 


- (g + a) 


? 

< 







-4a — (q — a) 2 


- (g + a) 


? 

< 






which is indeed the case. □ 
Lemma B.3 (details of the proof of Lemma |3.7|) . For a G N, q > 1 we have 

1 \ 2 4 



gi (g + l,a; 



) i afi) < ^T-(9 + l)((g + l) + « + l) 2 7£i 0) (B.3) 

52(9 + 1, «) \ (a,0) , 4 2 ( a ) 

Proof. We start with the bound ()B.3j) . We compute 

1 (a ,o) 1 = (2(g + 1) + a - l) 2 (2(g + 1) + a) 2 (2(g + 1) + a + 1) 

( 5 i (g + 1, a)) 27g (g + + 2 + a)^'? " ( 2 « + « + *)(% + X ) + 2 «) 2 (« + 2 + «) 2 (« + 1) 
1 (2g + l + a)(2(g + l)+a) 2 (2(g + l) + a + l) _ 4 (2g + 1 + a)(2g + a + 2) 2 (2g + a + 3) 



< 



q + 1 (2(g + l) + 2a) 2 (g + 2 + a) 2 g+1 (2g + 2a + 2) 2 (2g + 2a + 4) 2 

4 



_ + l 

We now turn to the bound (|B.4p . We first compute 

52(9 + 1,«) 2a (2(g + 1) + a - l)(2(g + 1) + a) (2g + a - l)(2g + a) 

51 (g + 1, a)gi(g, a) ~ (2(g + 1) + a - 2)(2(g + 1) + a) 2(g + 1) + 2a 2g + 2a 

_ 2a(2g + a + l)(2g + a-l) 
~ (2g + 2a + 2)(2g + 2a) 
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Then 



g 2 (q + l, a) Y%-i 



gi (q + l,a) gi (q,a)J J«fl) (q + l)(g + 2 + af 



lq+l 

2 



< 



< 



/2a(2g + a + l)(2g + a - 1) V 2(g + 1) + g + 1 1 

V (2g + 2a + 2)(2g + 2a) J 2(g - 1) + a + 1 (q + l)(q + 2 + a) 2 

4 (2a) 2 (2q + a + l) 2 (2q + a-l) 2 2q + a + 3 1 
g + 1 (2g + 2a + 2) 2 (2g + 2a) 2 2g + a - 1 (2g + 2a + 4) 2 

4 (2a) 2 
q + 1 (2g + 2a) 2 
4 



+ 1 



□ 

Lemma B.4. Zei £7 G C 1 (— 1, 1) and Zei (1 — x) a U(x) as well as (1 — x) a+1 U' be integrable. Further- 
more, let 

lim (1 - x) 1+a U(x) = and lim (1 + x)U(x) = 0. 

X— ¥l X— >— 1 

Consider hi, hi and /13 /rom 13. gj) . We define 



u q := J (1 - x) a ^(x)P (? (a ' 0) (x)(ix, 

6 3 := ^ (1 - x) a U\x)P { q afi) {x)dx. 

Then for q > 1 and a £ No i/ie following relationship holds: 

u q = hi(q,a)b q+ i + h 2 (q,a)b q + h 3 (q, a)6 g _i. 
Proof. From ([O]) we have for x — )• — 1 



and for x — > 1 



^ (1 - t) a p( a ' 0) {t)dt = 0(1 + x) 
j X (1 - t) a P^'°) {t)dt = 0{{\- x) a+l ) . 



Hence, using the stipulated behavior of U at the endpoints, the following integration by parts can be 
justified: 



J 1 (l-x) a U(x)P^(x)dx 
U{x) J X (l-t) a P^'°\x) 



- j 1 U'(x) j X (l - t) a p( a '°\t)dtdx. 



In particular, we note that b q is well-defined. Furthermore, 

1 u'(x) r {i-t) a p( a ' o) (t)dtdx 



1 



(1 - x) a U\x) (h 1 (q,a)P^(x) + h 2 (q,a)P^°\x) + h 3 (q,a)P^f\x)) dx 

= hi{q,a)b q+ i + h 2 (q,a)b q + h 3 (q,a)b q -i, 
where in the second equation we appealed to Lemma 13.21 (i). □ 
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Lemma B.5. For > -1 and U G C x (0, 1) n C((0, 1]) i/iere holds 

jf * x' 3 |C/(x)| 2 cfx < (-^-^ 2 jf 1 x^+ 2 |f/'(x)| 2 dx + ^^(l)! 2 - 
Proof. We define according to the notation of [TBI Thm. 330] 



/(*) :- 

For /3 > -1, HQ Thm. 330] states 



17' (x) 0<x<l , „, v /" 00 N , 10 x>l 

v ; and F x) := / / (t) dt = < 

x>l Jx \U(1)-U(x) 0<x<l. 



x^|F(x)| 2 dx< [-J-^Yj ^ +2 |f(x)| 2 dx 



Hence, 



1 x /3 (^7(x) - U{l)fdx < (-^—) C xP +2 \U'(x)\ 2 dx. 
By rearranging terms, we get 

-1 / o \ 2 /•! H 



x /3 |^7(x)| 2 dx < f-^-) / x /3+2 |?7 , (x)| 2 dx + |C/(l)| 2 / x^dx 

<^)Y* 2+VWI ^ + ^ |t,(1)|2 ' 



□ 



B.2 Selected proofs for Section [4] 

Lemma B.6 (details of Lemma l4.2j) . Let D be the Duffy transformation and T := T 2 x { — 1}. Then 
D(T) = r and D is an isometric isomorphism with respect to the L 2 (T)-norm. 

Proof. Obviously D is an isomorphism and by definition 77(r) = T, so we will only show the isometry 
property. 

Let u be a quadratic integrable function on T 3 and consider the transformed function u = u o D. 
We have 

1 r i 

2 



r j-u-i 
1 ' ' / (l + r ?1 )(l + r ?2 ) ^ 2 

« ( o 

1 J-l 



l^ 2 (r) = \u(m,V2,m)\ dr] 1 dr} 2 dr]3= J j \u(D(rn, n 2 , -1)) | dr/idr/ 2 

= J m^i,6,-i)i 2 rf66 = ii«ni2(r). 

□ 

Proof of Lemma {473[ With the definition of 77 _1 and the abbreviation n pq := 2p + 2q + 2 we get 

= pr ] (-2 ^| - 1) n (2p+i - o) ( 2 ^§ - 1) ^' o) fe) f 1 " 2 ^" 1 )' (^) P+9 • 
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Expanding P p (0 ' 0) (x - l)P 9 (2p+1 ' 0) (y - 1) = £Lo £? =0 <WV leads to 

-Z^Z^ CM ^(i + ^l) (1+6) , )k Pr 

fc=0 Z=0 VSi! W 

Since p^ npq ^ is a polynomial of degree r, we see by the last 2 terms in the sum above that ip Pl q jr G 

Pp+q+r(T 3 ). 

To see the orthogonality property, we transform to the cube S 3 and make use of (|3.ip three times 

r 3 



^p,qAO^p',q',r'(O d C = j ^ ^qAV^p'.q'^'iv) ( ~ ( ~ j^ J dr ? 

(i \ p+p'+l 

^) P^ +1 '% 2 )Pf' +1 '%2 



2 

2 



T V2- (2p+1) /' /V - ? 72) 2p+1 Pi 2p+1 ' 0) (ry2)Pf +1 '%2 



2p + 

2p+g+q'+2 



(-1 \ Zp+q+q +Z 

2 -V o Jo ±0 V^ l\i-m) npq P^ qfi \m)P^ pqfi \m)d m 



2p+l 2p + 2<7 + 2 

2 2 2 



2p+l pp 2p + 2q + 2 qq 2p + 2q + 2r + 3' 

□ 

Proof of Corollary \4-9\ To prove this corollary we want to make use of Lemma 13.31 Therefore, we 
have to clarify that the conditions in the lemma are satisfied. We proceed in two steps. First, we 
require u G C°°(T 3 ) and show the statement in this case and then we argue by density to achieve 
results in // 1 (T 3 ). 

Step 1: By assuming that u G C°°(T 3 ) we get u G C 1 (S 3 ). Hence, for fixed p and q, if we recall the 
definition of U Ptq , we see that the map 773 h-> U p AV3) 1S smooth on [—1, 1]. Considering the definition 
of U Pjq 

n (r, \ - u p>i^) 



{i-m) p+q 

we see that U Pt q G C 1 ([— 1, 1)) and that U Piq has at most one pole of maximal order p + q at the point 
773 = 1. In view of these preliminary considerations we conclude that the following limits exist and 
that the conditions in Lemma 13.31 are satisfied: 



and 



lim (1 - V3) 2p+2q+ '%AV3) = limfl - r, 3 ) P+q+3 U P Av 3 ) = 0. 

T73— >1 ' 



hm (1 + mPpAm) = 0. 
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Now the statement follows directly from Lemma 13.31 when looking at the definition of u Pt q jr an d u' pqr 
and consequently replacing U with U p>q and a with 2p + 2q + 1 . 

Step 2: Let u 6 ff 1 (7" 3 ). Since C^T" 3 ) is dense in H l (T^)-, there exists a sequence (u n ) ng N C 
C°°(T 3 ) such that it n — >• u in -ff^T" 3 ) for n — > oo. Because we have already proved that u n , n £ N 
satisfies our statement, ensuring that the sequences of coefficients u pqr and u'^ qr corresponding to 
u n converge for fixed p,q and r will conclude the proof: 

We have 

Since (ip P ,q,r)p,q,reNo forms an orthogonal basis for L 2 (T 3 ) and since i? 1 (T 3 ) C L 2 (7" 3 ), the maps 
F : u ^ Up t q jT are continuous linear functionals on H l (T^) and thus lim n _ > . 00 F(u n ) = F(u). 
In case of u' r we study the functionals F : u \-t u' Jhq r that map C°°(T 3 ) into R. Since F is a 
linear functional that is continuous with respect to the // 1 (7 _3 )-norm, we see by density of C°°(7" 3 ) 
in H l (T^) that it is indeed a well-defined continuous linear functional on H l (J~^) and thus again 
limn-^oo F(u n ) = F(u). □ 



C Extended versions of the tables of Section [2 
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N 


^^UeVoM 111, II n 111, II 1 


su PueV? N ii«ii 2 
ii i^i 


i 
i 


i.ioio4yioo04iyy 


U.o ( OUUUUUUUUUUU 


o 
z 


i.ozyoziizy /ybo / 


1 1 /19fi1 O091 fi771 8 

l.Moblzoolb / /lo 


Q 
O 


o 1 £0707£0900/i 1 ^ 

z.ioz (o i byoyu4io 


1 1 ^0700180^1 8r^o 
l.loll/ zU4ozblooz 


A 

4 


z.o4iUDzoy / loouy 


1 19C^90^000^/t^^7 

l.loooobUUob4ob / 





z.4oy4oyyizob4u / 


1 110007C8900917 

i.nyyz / ooooooi / 


1U 


O 701 (17fifiS8S0/l on 

z. / Ziy /bboooz4oU 


1 neofi70Q9^n7osfi 
l.Uozb / zoooU /yob 


10 


z.ozzlUoooUooooo 


l.UboOooolbUo4 10 


on 
zl) 


z.o / 4UD41ozzoyol 


1 o£ 1 1 1 o^ /i 7£ /i QQa 
!.UblllUb4 /b4oob 


o£ 
ZD 


z.yuoiz4oob4oiio 


1 0^^r^98919C010^ 

!.Uoboooo4ooU4yb 


on 
oil 


z.yzo4UolUzob4UU 


1 0^9190^/100001 O 

l.Uoo4oyo4zyUUlo 


OO 


z.yoy4oioUo4b / o4 


1 pit. 1 onnooQ7i /i o/i 
i.uoizuuyzo 1 i4y4 


/i n 
4U 


O 0/1071 1 OOOGKOQG 

z.y4y / iizyzybzoy 


i.U4youoUzoouiyz 


/I K 


z.yo / 4iioyb /uzui 


1 O/l Q1 Q00007K1 07 

i.u4oiozyyy /biz / 


p;o 
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O 0^9970707801 A O 

z.yboo /z /y / oyi4u 


1 0171 1 7^0C70^1^ 

1.1)4/ 11 /byo/yb4b 


00 


o ofisnoc;/! 7qoi i 

z.ybouyo4 /ouiio4 


1 O/l fiO/l 0C;70KQQ07 

l.U4bz4z0 /ybooz / 


£o 
bl) 


O 071 00000/171 1 £Q 

z.y / iyuyzU4 / iioo 


1 01^^100000^11 

l.l)4oo lUoyUoollo 




9 Q75D3Q9fi1 31 73f) 




70 


2.97764417211453 


1.04435662477469 


75 


2.97983833781387 


1.04389338318215 


80 


2.98170613800160 


1.04348732490767 


85 


2.98331100621813 


1.04312847760604 


90 


2.98470143645119 


1.04280906028005 


95 


2.98591505801301 


1.04252291272595 


100 


2.98698146107879 


1.04226509441461 


105 


2.98792419448666 


1.04203159687688 


110 


2.98876220322068 


1.04181913380901 


115 


2.98951087919110 


1.04162498545355 


120 


2.99018284042270 


1.04144688155757 



Table 3: Computed constants C for ID maximization problems 
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N 


II TT II '2. 

11 njvu| W) 


II TT „ . 1 1 ^ 

l|n^«|| £2(r) 


sup, 


1 1 TT ,», W'Z 

IPwu|| HlCT9) 


SUp u p-p II II || | 


SU P«GP 2 JV Hull* . , . 


Ll—T^o ^^ Mn,ll2 (AT I 1 \ 

«=/ 2 jv \\u\\ h1{t2) {n+i) 


1 


1.841709179979923 


1 


471718130438879 





630765682656827 


2 


2.482008261007026 


1 


705122181047698 





534605787746493 


3 


2.840187661660685 


1 


715716151553367 





502563892227422 


4 


3.069499343400879 


1 


698817794545474 





469837538241653 


5 


3.221409644442824 


1 


681456018543639 





451490548981791 


6 


3.328227344860761 


1 


668376974400849 





442841777476419 


7 


3.407909278234681 


1 


658518924126468 





439947515744027 


8 


3.470131098161663 


1 


650883702558015 





438357597812212 


9 


3.520381327094282 


1 


644802661753009 





437324308925344 


10 


3.561974205023130 


1 


639847070385123 





436649057543612 


11 


3.597047215946491 


1 


635732061317178 





436162962404347 


12 


3.627055403122174 


1 


632261308353419 





435795222777319 


13 


3.653032616328430 


1 


629295122853766 





435497902586582 


14 


3.675739731999510 


1 


626731407637615 





435247417526015 


15 


3.695751965776597 


1 


624493814472113 





435028483384906 


16 


3.713514379154265 


1 


622524114095057 





434833152302013 


17 


3.729377544445405 


1 


620777127516557 





434656212926501 


18 


3.743622134346961 


1 


619217268307188 





434494915639700 


19 


3.756475749181842 


1 


617816128453177 





434347360750190 


20 


3.768125200720072 


1 


616550757911536 





434212531866692 


21 


3.778725310656065 


1 


615402415968316 






22 


3.788405606490951 


1 


614355650089980 






23 


3.797275282737163 


1 


613397606293759 






24 


3.805427098445950 


1 


612517505912479 






25 


3.812940340675182 


1 


611706243752921 





433709285223336 


26 


3.819883212382239 


1 


610956076034725 






27 


3.826314685175567 


1 


610260375564773 






28 


3.832286023912474 


1 


609613437845351 






29 


3.837841988875210 


1 


609010326178230 






30 


3.843021842931179 


1 


608446746919096 





433462666731238 


31 


3.847860156122867 


1 


607918948256749 






32 


3.852387490412042 


1 


607423637503739 






33 


3.856630952540156 


1 


606957913069210 






34 


3.860614671188471 


1 


606519208163681 






35 


3.864360185593495 


1 


606105243942980 





433441925213746 


36 


3.867886785255372 


1 


605713990297966 






37 


3.871211788573052 


1 


605343632872021 






38 


3.874350789295843 


1 


604992545187163 






39 


3.877317859826744 


1 


604659264974558 






40 


3.880125733031547 


1 


604342473993514 





433618345964902 


41 


3.882785952906251 


1 


604040980750689 






42 


3.885309010715330 


1 


603753705645651 






43 


3.887704458222425 


1 


603479668157049 






44 


3.889981011011230 


1 


603217975745874 






45 


3.892146634658522 


1 


602967814218030 





433964716705375 


46 


3.894208624113123 


1 


602728439323389 






47 


3.896173670029407 


1 


602499169413209 






48 


3.898047920431772 


1 


602279379001790 






49 


3.899837032314585 


1 


602068493106013 






50 


3.901546220032895 


1 


601865982255488 





434455030379973 


51 


3.903180295826488 


1 


601671358082506 






52 


3.904743708152147 


1 


601484169413596 






53 


3.906240573788223 


1 


601303998798579 






54 


3.907674708463756 


1 


601130459422183 






55 


3.909049652499364 


1 


600963192345851 





435064216731618 



Table 4: Extended version of Table [2] 
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